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The holographic principle in a radiation dominated universe, as discussed 
first by Verlinde [hep-th/008140], is extended so as to incorporate the case of 
a bulk-viscous cosmic fluid. This corresponds to a non-conformally invariant 
theory. Generalization of the Cardy-Verlinde entropy formula to the viscous 
case appears to be formally possible, although on physical grounds one may 
' question some elements in this type of theory, especially the manner in which 

the Casimir energy is evaluated. Also, we consider the observation made by 
Youm [Phys. Lett. B531, 276 (2002)], namely that the entropy of the universe 
is no longer expressible in the conventional Cardy-Verlinde form if one relaxes 
the radiation dominance equation of state for the fluid and instead merely 
. assumes that the pressure is proportional to the energy density. We show that 

""q ' Youm's generalized entropy formula remains valid when the cosmic fluid is no 

Q" 1 longer ideal, but endowed with a constant bulk viscosity. In the introductory 

part of this article, we take a rather general point of view and survey the 
bJT), essence of cosmological theory applied to a fluid containing both a constant 

shear viscosity and a constant bulk viscosity. 
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1 Introduction 



Perfect fluid models have for a long time been used in cosmological theory. The 
introduction of the viscosity concept came later, and there has in fact been an open 
question as to whether the viscosity concept is needed for an explanation of the 
observed quantities in the universe. Misner was probably the first to introduce 
viscosity in cosmology in connection with his study of how initial anisotropics in 
the early universe became relaxed. Cosmological model with viscosity have later 
now and then been discussed in the literature, from various points of view. A useful 
review of the subject, with many references to the earlier literature up to 1990, was 
given by Gr0n |2j. 
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From a physical point of view it would in our opinion be almost surprising if the 
viscosity concept were not of importance in cosmology. An essential ingredient of 
the Friedmannian model of the universe is after all to borrow the energy-momentum 
tensor from non- viscous fluid dynamics and insert it into the Einstein equations. 
From fluid dynamics we know that in many situations the non-viscous theory is 
inadequate, in particular, if anisotropies or turbulence effects are involved. Use of 
the molecular shear and bulk viscosities means in effect an expansion of the theory 
to first order in the deviations from thermal equilibrium. Moreover, once turbulence 
occurs in the velocity field it is the Reynolds stresses, rather than the stresses caused 
by molecular viscosities, that govern the behavior of the fluid. In view of these 
known facts from ordinary fluid dynamics it is natural to be somewhat reluctant to 
regard cosmology as an exceptional case for which the viscosity concepts are of no 
use. One main reason why the viscous cosmological theory has so far not gained 
large attention is undoubtedly that the viscosity-generated entropy a per baryon is 
observed to be very large, a ~ 10 9 , and previous investigations have shown that a 
straightforward use of the kinematically derived bulk viscosity in the isotropic and 
homogeneous universe is unable to explain this large magnitude EJ IS] . 

The fluidlike behavior of the universe during its early epochs is in all probability 
very complicated. One particular facet of the problem is the existence of one or more 
phase transitions that may have led to sudden changes whose description is naturally 
given in terms of the phenomenological viscosity coefficients. An eclatant example of 
this sort is the transition from the de Sitter universe back to the FRW universe at the 
end of the inflationary era, at t ~ t~ 33 s. We may call the effective viscosity operative 
during such a phase transition (bulk viscosity in case of an isotropic universe) an 
"impulsive" viscosity. Examples of this kind of theory were examined in Refs. [Sill]- 

Another interesting aspect of viscous cosmology that has arisen recently is the 
influence from shear viscosity. The shear viscosity comes into play in connection 
with anisotropy. Thus Weinberg |5] derived an analytic formula for the traceless 
part of the anisotropic stress tensor due to freely streaming neutrinos. Another 
recent example is the suggestion of Kovtun et al. [5] and Karch [101 (cf. also the 
discussion in that there exists in cosmology a universal lower bound on the 
ratio rj/s, where n is the shear viscosity and s the entropy per unit volume. 

The main topic to be dealt with in the present article is the incorporation of 
the bulk viscosity £ in the so-called holographic principle for the early universe. 
In the case of a radiation dominant ideal-fluid universe, Verlinde ^2] put forward 
the idea that there exists a bound on the subextensive entropy associated with the 
Casimir energy. This bound is called the holographic bound. When this bound is 
saturated, there exists a formal coincidence between the Friedmann equation for 
H 2 = (a/ a) 2 and the Cardy entropy formula known from conformal field theory 
|13| I14j . The question that naturally arises is whether this merging between the 
holographic principle, the entropy formula from conformal field theory, and the 
Friedmann equation from cosmology, is of deeper physical significance and thus 
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reproducible under more general conditions, or if it is just a formal coincidence. 
As one would expect, the Verlinde proposal has been the subject of study in cases 
where more general effects have been accounted for. For instance, Wang et al. |15j 
have considered universes having a cosmological constant, and Nojiri et al. |16j 
have considered quantum bounds for the Cardy- Verlinde formula. There are several 
papers in this area of Padmanabhan et al. ^7j, and of Cai et al. the paper of 
Youm contains an extensive list of references. We also mention the considerable 
interest that has arisen in connection with entropy and energy as following from 
quantum and thermal fluctuations in conformal field theories |20| . 

We begin in the next section by giving a survey of cosmological theory applied 
to a fluid that is general enough to possess both a constant shear viscosity r\ and 
a constant bulk viscosity £. Thereafter, in Sect. 3 we consider the Verlinde setting 
when the fluid is taken to possess a bulk viscosity only. In Sect. 4 we carry out 
an analogous reasoning for the case that the bulk-viscosity fluid is not necessarily 
radiation dominated but instead required to obey the equation of state 



with 7 a constant. Sections 3 and 4 are based upon material given previously in 
Refs. (23 [22l |2H]. We supply the analysis with some numerical estimates; thus 
we exploit the smallness of £ in the plasma era to calculate the viscosity-induced 
correction to the scale factor a(t), assuming for simplicity that the spatial curvature 
k is zero. In Sect. 5 we test the proposed holographic entropy bound on an example, 
taken from the plasma era in the early universe. 

We mostly use natural units, with h = c = kg = 1. 

2 Survey of viscous cosmology 

We use the convention in which the Minkowski metric is (-++-1-). Let = (U°, U l ) 
be the four-velocity of the cosmic fluid. In comoving coordinates, U° = 1, U l = 0. 
Let g^u be the general metric. Using the projection tensor 



P = (7 - 1)P, 



(1) 



h 



(2) 



we define the rotation tensor as 




(3) 



and the expansion tensor as 



0, 



\{u^K + u v;oi h«). 



(4) 



The scalar expansion is 9 = Oft = U^- 



The shear tensor, as defined by 



0", 



{IV u fJ,U 




(5) 
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is traceless, <r^ = 0. The following decomposition of the covariant derivative is often 
useful: 

1 

U^=Ld^ + a flI/ + -h fl „e-A^U„, (6) 

where An = U„ = U V U^- U is the four-acceleration of the fluid. 

We write down the expression for the energy-momentum tensor T^ u of the vis- 
cous fluid, taking into account also the conduction of heat. As already mentioned, 
rj and £ are respectively the shear and the bulk viscosities; if moreover k is the 
thermal conductivity (all quantities taken in accordance with their nonrelativistic 
definitions), then 

Q» = - K h^(T„ + TA v ) (7) 
is the spacelike heat flux density four-vector, and 

= P U^U U + (p- (9) V " 2r/<v + Q„U V + Q U U^. (8) 

Here p is the mass-energy density and p is the isotropic pressure, both taken in the 
local rest inertial frame. The last term in Eq. (|7|. containing TA V , is of relativistic 
origin. If one ignores this term, one is left with = —kWT u . This expression 
is defined such that in a local rest inertial frame (designated with a "hat") Qq = 
0, whereas Q- = — acTj is the heat energy per unit time crossing a unit surface 
orthogonal to the unit vector ej. 

Consider next the production of entropy. It is here instructive to start from 
nonrelativistic theory. If itj are the nonrelativistic velocity components, and a the 
nonrelativistic entropy per particle (baryon), then the ordinary entropy per unit 
volume is S = na, n being the baryoun number density. From nonrelativistic theory 
we have 

f = ^(*« " \W ■ u) 2 + i (V • u) 2 + ^(VT) 2 , (9) 

where = Uu k y The transition to relativistic theory may be made via the effective 
substitutions 

Oik -> 6^u, 5 ik -> h^, V • u -> 9, -nT k -> (10) 
from which we obtain 

S% = ^<j^ v + ^ 2 + Q^. (11) 
Here, is the entropy four-vector 

= noU* + -Q*. (12) 

Our treatment above follows Ref . [0] • The same result is obtained from a more careful 
analysis taking into account the relativistic thermodynamic equations H3] . 
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In the case of thermal equilibrium, = 0. Moreover, in accordance with 
usual practice we omit the shear viscosity in view of the assumed complete isotropy 
of the cosmic fluid, although we have to remark that this is actually a nontrivial 
point. The reason is that the shear viscosity is usually so much greater than the bulk 
viscosity. Typically, after termination of the plasma era at the time of recombination 
(T ~ 4000 K), the ratio n/C, as calculated from kinetic theory is as large as about 
10 12 0. Thus, even a slight anisotropy in the fluid would easily outweigh the effect 
of the effect of the very small bulk viscosity. 

3 Cardy-Verlinde formula in a bulk-viscous radiation 
dominated universe 

Assume now that the metric is of the Friedmann-Robertson- Walker (FRW) type, 
ds 2 = _ dt 2 + a 2 (t) + r 2 dQ^j , (13) 

k = —1,0, 1 being the curvature parameter. In comoving coordinates, U^ v = T^,, 
so that in view of the standard relations for the Christoffel symbols we get U^- u = 
h^a/a, with d = da/dt. From these equations we see that the rotation and shear 
tensors both vanish, 

w M „ = °> = °> ( 14 ) 

whereas the scalar expansion is 

= 3a/a = 3H, (15) 

H being the Hubble parameter. The four-acceleration is zero, = 0. The energy- 
momentum tensor becomes 

Too = p, T ok = 0, T ik = (p - Q0)g ik . (16) 

There is thus no conduction of heat in the FRW space, which is a homogeneous 
space. The entire effect of the bulk viscosity is to reduce the pressure by an amount 
£8, so that the effective pressure becomes p = p — (9. 

When applied to the FRW space, Eqs. ((Til and JEJ) yield 



S% = ^# 2 , S° = na, S* = 0. (17) 



Consider now the Einstein equations 



Rfj,u ~ \.Rg^v + ^9nv = 8vrGT^, (18) 
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from which we obtain the first Friedmann equation (the "initial value equation") 

tt2 8itG A k 
3 3 a z 

p = E/V being the energy density. This equation contains no viscous term. The sec- 
ond Friedmann equation (the "dynamic equation"), when combined with Eq. ([190. 
yields 

H = _4vrG(p + p) + 4, (20) 
or 

in which the presence of viscosity is explicit. 

The differential conservation for energy, T 0u ]U = 0, yields 

P+(p + P )9 = (6 2 . (21) 

The conservation equation for baryon particle number is 

(nU^ = 0, (22) 

which means that na 3 = constant in the comoving frame. Then from Eq. (|17j) 
we obtain, when substituting from Eq. ([12|) and observing that Q° = in the 
comoving frame, 

no = ^6 2 . (23) 

We now recall that the entropy of a (1+1) dimensional CFT is given by the 
Cardy formula [131 114j 




S = 2jtJ- [Lq-— , (24) 



where c is the central charge and Lq the lowest Virasoro generator. 

Let us assume that the universe is closed, and has a vanishing cosmological con- 
stant, k = +1, A = 0. This is the case considered by Verlinde 12 (in his formalism 
the number n of space dimensions is set equal to 3). The Friedmann equation l[19j) 
is seen to agree with the CFT equation (|24|) if we perform the substitutions 

L °-l E «- s -fr < 25 > 

These substitutions are the same as in Ref.|12j. We see thus that Verlinde's argument 
remains valid, even if the fluid possesses a bulk viscosity. Note that no assumptions 
have so far been made about the equation of state for the fluid. At this point the 
following question thus naturally arises: is there a deeper connection between the 
laws of general relativity and those of quantum field theory? 

Continuing this kind of reasoning, let us consider the three actual entropy defi- 
nitions. 
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• First, there is the Bekenstein entropy |26j . 

2-TT 

S B = —Ea. (26) 

The arguments for deriving this expression seem to be of a general nature; in 
accordance with Verlinde we find it likely that the Bekenstein bound S < Sb 
is universal. We shall accept this expression for Sb in the following, even when 
the fluid is viscous. 

• The next kind of entropy is the Bekenstein-Hawking expression Sbh, which is 
supposed to hold for systems with limited self-gravity: 

Sb » = m- (27> 

Again, this expression relies upon the viscous-insensitive member ()19|) of Fried- 
mann's equations. Namely, when A = this equation yields 

Sb < Sbh when Ha < 1, 

S B > Sbh when Ha > 1. (28) 

The borderline case between a weakly and a strongly gravitating system is 
thus at Ha = 1. It is reasonable to identify Sbh with the holographic entropy 
of a black hole with the size of the universe. 

• The third entropy concept is the Hubble entropy Sh- It can be introduced by 
starting from the conventional formula A/ AG for the entropy of a black hole. 
The horizon area A is approximately H~ 2 , so that 

IT — 2 ffV 

since V ~ H~ 3 . Arguments have been given by Easther and Lowe |27j . 
Veneziano |28j . Bak and Rey and Kaloper and Linde [SO] for assum- 
ing the maximum entropy inside the universe to be produced by black holes 
of the size of the Hubble radius (cf. also According to Verlinde the FSB 

prescription (see ^2] f° r a closer discussion) one can determine the prefactor: 

SB = %■ (30) 

It is seen to agree with Eq. (|25j). 

One may now choose (see below) to define the Casimir energy Ec as the violation 
of the Euler identity: 

E c = 3(E + pV- TS) (31) 
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where, from scaling, the total energy E can be decomposed as {Ee is the extensive 
part) E(S, V) = E E {S, V) + \E C {S, V). Due to conformal invariance the products 
Ee a and Eq a are independent of the volume V, and a function of the entropy S 
only. From the known extensive behaviour of Ee and the sub-extensive behaviour 
of Eq one may write (for CFT) 

E E = -^5 4 / 3 , E C = /-S 2 /3 ; (32) 

where a, j3 are constants whose product for CFTs is known: yja(3 = n = 3 (this 
follows from the AdS/CFT correspondence, cf. [12 ). From these expressions it 
follows that 

S = ^^E C {2E-E C ). (33) 

This is the Cardy-Verlinde formula for the radiation dominated universe. Identifying 
Ea with Lq and Eq a with c/12 we see that Eq. (|33|) becomes the same as Eq. (|24[l. 
except from a numerical prefactor which is related to our assumption about n = 3 
space dimensions instead of n = 1 as assumed in the Cardy formula. 

The question is now: can the above line of arguments be carried over to the 
case of a viscous fluid? The most delicate point here appears to be the assumed 
pure entropy dependence of the product Ea. As we mentioned above, this property 
was derived from conformal invariance, a property that is absent in the case under 
discussion. To examine whether the property still holds when the fluid is viscous 
(and conformal invariance is lost), we can start from the Friedmann equations (|19|) 
and (|20|> . in the case k = 1, A = 0, and derive the "energy equation", which can be 
transformed to 

±( P a 4 ) = ( P -3p)a 3 . (34) 
Thus, for a radiation dominated universe, p = p/3, it follows that 

|(pa 4 ) = C^V. (35) 

Let us compare this expression, which is essentially the time derivative of the volume 
density of the quantity Ea under discussion, with our earlier expression (|23|): both 
the two time derivatives are seen to be proportional to Since £ is small, we 
can insert for a = a{t) the expression pertinent for a non-viscous, closed universe: 
a{t) = a* sin 7/, where r\ here denotes conformal time, and a* is the constant 



(8irG/3)p in af n . (36) 

The subscript "in" designates the initial instant of the onset of viscosity. Imagine 
now that Eqs. (|35|) and ()23|) are integrated with respect to time. Then, since the 
densities C~ 1 pa 4 and (~ 1 na can be drawn as functions of t, it follows that pa 4 can 
be considered as a function of na, or, equivalently, that Ea can be considered as a 
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function of S. We conclude that this property, previously derived on the basis of 
CFT, really appears to carry over to the viscous case. 

The following point ought to be commented upon. The specific entropy a in 
Eq. (|23f) is the usual thermodynamic entropy per particle. The identification of 
S with HV/(2G), as made in Eq. (|25jl. is however something different, since it is 
derived from a comparison with the Cardy formula (|24jl. Since this entropy is the 
same as the Hubble entropy Sh we can write the equation as nan = H/(2G), where 
on is the Hubble entropy per particle. This quantity is different from <r, since it 
does not follow from thermodynamics plus Friedmann equations alone, but from the 
holographic principle. The situation is actually not peculiar to viscous cosmology. 
It occurs if £ = also. The latter case is easy to analyze analytically, if we focus 
attention on the case t — ► 0. Then, for any value of k, we have a oc i 1//2 , implying 
that H = l/(2t). Moreover, the equation of continuity {nU^)-a = implies, as we 
have seen, that na 3 = constant for a FRW universe, so that n oc i~ 3 / 2 . The above 
equation for an then yields an oc t 1 / 2 . This is obviously different from the result 
for the thermodynamic entropy a: from Eq. (|2M|) we simply get a = constant when 
£ = 0. The two specific entropies are thus different even in this case. 

3.1 Remarks on the physical interpretation 

Let us make three remarks on the interpretation of the above formalism. They are 
based on physical, rather than mathematical, considerations. 

• First, one may wonder about the legitimacy of defining the Casimir energy such 
as in Eq. Usually, within the Green function approach, in a spherical 
geometry the Casimir energy is calculated indirectly, by integration of the 
Casimir surface force density / = — (l/4:7ra 2 )dE/da. The force / in turn 
is calculated by first subtracting off the volume-dependent parts of the two 
scalar Green functions; this is in agreement with the physical requirement 
that f — > at r — ► oo. (The typical example of this configuration is that of 
a conducting shell; cf., for instance, Ref. |S2]-) That this kind of procedure 
should lead to the same result as Eq. Q31|). which merely expresses a violation 
of the thermodynamic Euler identity, is in our opinion not evident. 

• Our second remark is about the physical meaning of taking the Casimir energy 
Ec to be positive. Verlinde assumes that Eq is bounded by the total energy 
E: Eq < E. This may be a realistic bound for some of the CFTs. However, 
in general cases, it is not true. For a realistic dielectric material it is known 
that the full Casimir energy is not positive; the dominant terms in Ec are 
definitely negative. From a statistical mechanical point of view this follows 
immediately from the fact that the Casimir force is the integrated effect of 
the attractive van der Waals force between the molecules. Now the case of a 
singular conducting shell is complicated - there are two limits involved, namely 
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the infinitesimal thickness of the shell and also the infinite conductivity (or 
infinite permittivity) - and a microscopical treatment of such a configuration 
has to our knowledge not been given. What is known, is the microscopical 
theory for a dielectric ball. Let us write down, for illustration, the expression 
derived by Barton [3S] for a dilute ball: 

2 ( 3 V 7 A 11 23 1 



+ 7 { 128tt 2 A 4 + 360vr3 A3 20vr 2 A + 1536vr a) ' 

where 7 = (e — l)/e, A is the surface area, and A is a cutoff parameter. This 
expression, derived from quantum mechanical perturbation theory, agrees with 
the statistical mechanical calculation in Ref. |34| . and also essentially with 
Ref. |35j (there are some numerical factors different in the cutoff-dependent 
terms). It is evident from this expression that the dominant, cutoff, dependent 
voume terms, are negative. 

We see that there remains one single, cutoff independent, term in Eq. IJ37JI . 
This term is in fact positive. It can be derived from macroscopic electrody- 
namics also, by using either dimensional continuation or zeta-function regu- 
larization, as has been done in Ref. [SHI- In the present context the following 
question becomes however natural: how can a positive, small, cutoff depen- 
dent term in the Casimir energy play a major role in cosmology? In another 
words, why should the matter necessarily be conformal? Of course, our uni- 
verse is different from a dielectric ball, and we are not simply stating that 
Verlinde's method is incorrect. Our aim is merely to stress the need of cau- 
tion, when results from one field in physics are applied to another field. In any 
case, all this suggests that the consideration of non-conformally invariant sit- 
uations should significally change the dynamical entropy bounds and bounds 
for Casimir energy. 

Our discussion on the generalized Verlinde formula in the present section was 
based upon the set of cosmological assumptions {p = p/3, k = +1, A = 0}. 
The recent development of Wang et al. ^H] is interesting, since it allows for 
a nonvanishing cosmological constant (still assuming a closed model). One of 
the scenarios treated in [15] is that of a de Sitter universe (A > 0) occupied 
by a universe-sized black hole. A black hole in de Sitter space has the metric 

ds 2 = -f{r)dt 2 + f^^dr 2 + r 2 dQ 2 , (38) 

where f(r) = 1 — 2MG/r — Ar 2 /3. The region of physical interest is that lying 
between the inner black hole horizon and the outer cosmological horizon, the 
latter being determined by the magnitude of A. 
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Although we do not enter into any detail about this theory, we make the 
following observations: the above metric is static; there is no time-dependent 
scale factor involved, and the influence from viscosity will not turn up in the 
line element. Moreover, Wang et al. make use of only the member (|19|) of 
Friedmann's equations which, as we have noticed, is formally independent of 
viscosity. 

Does this imply that viscosity is without any importance for the present kind of 
theory? The answer in our opininon is no, since the theory operates implicitly 
with the concept of the maximum scale factor a max in the closed Friedmann 
universe. In order to calculate a max , one has to solve the Friedmann equation 
(|2U|) also, which contains the viscosity through the modified pressure p. Thus, 
viscosity comes into play after all, though in an indirect way. 



4 Cardy-Verlinde formula in the presence of a general 
equation of state 

The paper of Youm mentioned above |19j is interesting, since it shows that the 
entropy of the universe can no longer be expressed in the conventional form of 
Eq. 1)33(1 if one relaxes the radiation dominance state equation and instead assumes 
the more general equation And this brings us naturally to the problem of 

how the presence of a bulk viscosity, together with Eq. influences the entropy 
formula. Actually, it will turn out that the modified entropy formula found by Youm 
still persists, even when a constant bulk viscosity is allowed for. The central formula 
is Eq. (g21 ) below. 

The main part of the formalism has been given already. Thus Friedmann's 
equations take the same form (|19|) and 1|2U[) as before (with k = +1,A = 0), the 
conservation equation for energy is still Eq. (|21j) . and also the production rate for 
entropy is as in Eq. (|23[). However, Eq. (|3*5|) gets replaced by 

j t {pa^) = C9 2 a^. (39) 

We can now carry out the same kind of reasoning as above: Since £ is assumed 
small, we can use for a = a(t) the same expression as for a nonviscous closed universe, 
namely a{t) = a* sin 77 with a* given by Eq. (|3l))) (although this approximate expres- 
sion strictly speaking assumes that the universe is radiation dominated). Imagine 
that Eqs. (j3H|) and (|23j) are integrated with respect to time. Since C~V« 37 and 
Q^na can be drawn as functions of t, it follows that pa 37 can be considered as a 
function of na. Then, since the total energy is E ~ pa? and the total entropy is 

5 ~ naa 3 , it follows that Ea 3 ^~^ is independent of the volume V and is a function 
of S only. This generalizes the pure entropy dependence of the product Ea, found 
by Verlinde ^2] hi the case of a nonviscous radiation dominated universe. And it is 
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noteworthy that the derived property of Ea 3 ^^ 1 ^ formally agrees exactly with the 
property found by Youm ^5] when £ = 0. 

Let us carry out the analysis a bit further, and write the total energy E as a 
sum of an extensive part Ee and a subextensive part Ec, as in the previous section. 
Under a scale transformation S — > XS and V — * XV with constant A, Ee scales 
linearly with A. But the term Eq scales with a power of A that is less than one: as 
Ec is the volume integral over a local energy density expressed in the metric and 
its derivatives, each of which scales as A -1 / 3 , and as the derivatives occur in pairs, 
the power in A has to be 1-2/3= 1/3. Thus we have 

E E (XS, XV) = XE E (S, V), E C (XS, XV) = A 1 / 3 E C (S, V), (40) 

which implies 

Together with the decomposition of E(S, V), this leads to 



S 



Ec(2E — Ec) 



3 

37-1 



(42) 



This is the generalized Cardy-Verlinde formula, in agreement with Eq. (20) in 
Youm's paper, reducing to the standard formula (with square root) in the case 
of a radiation dominated universe. In conclusion, we have extended the basis of 
Eq. (|42[) so as to include the presence of a constant bulk viscosity in the cosmic 
fluid. 



4.1 Numerical estimates for a spatially flat universe 

It is of physical interest to supplement the above considerations with some simple 
numerical estimates, showing, in particular, the order of magnitude of the viscosity 
terms. Let us assume, as mentioned above, that the effect of the bulk viscosity 
effectively sets in at some initial instant t = t; n and that £ is thereafter constant for 
the times that we consider. For definiteness we choose 

t in = 1000 s (43) 

after the big bang. The universe is then in the plasma (or radiation) era; it consists 
of ionized H and He in equilibrium with radiation. The particle density is n; n ~ 
10 19 cm -3 , and the temperature is T m ~ 4 x 10 8 K. The advantage of considering 
this relatively late stage of the universe's history is that the magnitude of £ can be 
calculated using conventional kinetic theory. At t = t{ n one finds 

C = 7.0 x 10~ 3 gem" 1 s^ 1 (44) 



(cf. |)] and further references therein). 
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For our estimate purposes it appears reasonable to assume that the influence 
from the spatial curvature is not very important. For simplicity let us put k = 0. 
This implies that there exists the following simple differential equation for the scalar 
expansion 6 : 



Ht) + \l8 2 {t) ~ 127rGC#(t) = 0, 
which can be solved to give the expression for the scale factor: 

2/(3 7 ) 



a(t) 



air 



1 



where 



t c = (127TGC) -1 . 

The corresponding expression when viscosity is absent, is 



a(i,C = 0) 



1 + ^M* - ti 



1 2/(37) 



(45) 

(46) 
(47) 

(48) 



The ratio between the expressions (|46|) and (|48j) reduces to unity in the limit when 
(t — ti n )/t c <C 1. This is we would expect. The influence from viscosity generally 
turns up only in the factor t c , and the effect becomes strengthened when t c becomes 
smaller, i.e., when £ becomes larger. The numerical values given above for the 
instant ii n = 1000 s correspond to 



1.5 x 10~ 3 s -1 , 



tr 



5.1 x 10 28 s. 



(49) 



4.2 Perturbative expansion for a radiation dominated closed uni- 
verse 

The smallness of £ makes it natural, as an alternative to the approach of the previous 
subsection, to make a Stokes expansion in £. For simplicity we now assume that 
7 = 4/3, i.e., that the universe is radiation dominated. We put k = 1. Let subscript 
zero refer to quantities in the nonviscous case. We write the first order expansions 
as 

a = a (l + C«i), P = Po(l + CPi), (50) 

where the functions a\ and p\ are of zeroth order in they are regarded as functions 
of t or alternatively as functions of the conformal time n. Correspondingly, the scalar 
expansion 9 = 3a /a is written as 

e = e (i + (e 1 ). (51) 

The Friedmann equation ()19|) . and Eq. ()39|) . now yield to the zeroth order 



2 o = 24irG Po --, 



(52) 
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Pinafn = Poao, (53) 

and to the first order 

e 2 e 1 = 12nGp oPl + ^, (54) 
a 

p in af n (pi + 4di) = 0q Oq. (55) 
The solutions of Eqs. (|52|) and (|53|) are 

3 1 1 
a = a* sin 77, p = — — — — , 
»7rG a% sin 77 



3 cos 77 

#0 = r-o— j 

a* sm 77 

where a* as before is given by Eq. ()36|) . From Eq. (|55|) we now have 



(56) 



dpi 4dai . 2 

— 1 : — = 247rGa* sm 77 cos 77, (57) 

an an 

which can be integrated from 77 = 7/i n onwards: 

Pi + 4ai = 87rGa^(cos 3 77^ — cos 3 77). (58) 
We have here assumed that a± = p\ = at the initial instant t = t- m . With 

1 = ^itan77 (59) 

GE77 

(cf. Eq. (|5U[)). we obtain from Eq. (|54|) a first order differential equation for the 
scale factor correction: 

sin 277 -^i + 2(1 + cos 2 77)01 + 8nGa* cos 3 77 = 87rGa* cos 3 77i n . (60) 
an 

After some calculation we find the following solution, again observing the initial 
conditions at t = t{ n : 



, , 4-irGa* 

ai{V) = . 2 
sm 7/ 



/ 1 2 1 

COS 7/in + ( - COS 2f] — COS 77; n — - COS 277; n ) COS 77 



(61) 



Once the scale factor correction a\ is known, the corresponding density correction 
pi follows from Eq. (|58|) . 

From Eq. ()61|) it is apparent that the relative correction £ai to the scale factor is 
of order AirGC,a* or, in dimensional units, AirGCa^/c 3 . We here note that AttGC/c 2 = 
l/(3i c ), thus about 6.5 x 10 -30 s _1 according to Eq. (|47|) . which in dimensional form 
reads t^ 1 = 12ttG(/c 2 . The quantity a*, according to Eq. (f5^|) . is the maximum 
value of the scale factor in a nonviscous k = 1 theory. Let us put a* equal to the 
commonly accepted value of the radius of the universe, i.e., a* = 10 28 cm. Then we 
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obtain 47rG£a*/c 3 = 2 x 1CP 12 . The relative correction to the scale factor is thus 
in this case very small. Physically, this is due to the fact that we are considering 
an instant relatively late in the history of the universe. If the bulk viscosity £ were 
higher at earlier times, or, if there were an anisotropic stage present in the early 
universe at which the enormously higher shear viscosity would come into play 
then the effect would be significantly enhanced. 



5 Remarks on shear viscosity and the holographic en- 
tropy bound 

We round off this article by making some remarks on the influence from shear viscos- 
ity in cosmology. The shear viscosity concept, as mentioned in the Introduction, has 
received considerable attention recently [HI El I10| 111] . What we shall focus attention 
on here, is the suggestion of Kovtun et al. that there exists a cosmological uni- 
versal lower bound on the ratio r//s, rj denoting as usual the shear viscosity and s the 
entropy per unit volume. The mentioned authors are concerned with the infrared 
properties of theories whose gravity duals contain a black brane with a nonvanish- 
ing Hawking temperature, the point being that the infrared behavior is governed by 
hydrodynamical laws. If we for definiteness consider a stack of iV non-extremal D3 
branes in type IIB supergravity, the metric near the horizon is given by 

ds 2 = ^[-f{r)dt 2 + dx 2 + dy 2 + dz 2 \ + -^jr-Ar 2 + R 2 dQ 2 , (62) 
R z r z j{r) 

where R oc iV 1 / 4 is a constant, and f(r) = 1 — r 4 /r 4 with ro being the horizon. The 
Hawking temperature of this metric is T = ro/irR 2 , and r] and s are given by 

rj = -ttN 2 T 3 , s = -tt 2 N 2 T 3 . (63) 
8 2 



Thus, when using hereafter dimensional notation, 

6.08 x KT ia Ks. (64) 



'/ h in -13 



s A-nks 

The conjecture of Kovtun et al. is that the value in Eq. (|64j) is a lower bound for rj/s. 
Since this bound does not involve the speed of light, the authors even conjecture 
that this bound exists for all systems, including non-relativistic ones. The idea has 
recently been further elaborated in ^0], arguing that the bound follows from the 
generalized covariant entropy bound. 

Let us check the proposed bound, by considering an example taken from the 
plasma era in the early universe. We choose the same instant as in subsection 
4.1, namely ti n = 1000 s, corresponding to n ~ 10 19 cm~ 3 (subscript "in" here 
omitted), T ~ 4 x 10 8 K, whereas the energy density is pc 2 = a T T^, where a r = 
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7r 2 kg I (157l 3 c 3 ) = 7.56 x 1CP 15 erg cm -3 K -4 is the radiation constant. The viscosity 
coefficients - whose existence is due to the imperfectness of thermal equilibrium - 
can be calculated from relativistic kinetic theory. Let x = m e c 2 /ksT be the ratio 
between electron rest mass and thermal energy; when x 3> 1 it is convenient to use 
the polynomial approximations |37| (cf. also (HJ) for the evaluation of rj and £: 

_ 5m6 c 8 C(3) _ 4 _ irc 2 h 3 n 3 
V ~ 9^h 3 e^n X ' C - 256e 4 C(3) X ' (b5) 

C(3) = 1.202 being the Riemann zeta function. At T = 4 x 10 8 K one has x = 14.8, 
leading to 

rj = 2.8 x 10 14 gcm^s" 1 , (66) 

which is enormously greater than £ given by Eq. (|44[). We note that both rj and £ 
contain /i. 

The entropy density, in view of the radiation dominance, is given by 

s = \a r T z = 6.45 x 10 11 erg cm" 3 K~\ (67) 

and so 

-= 435 Ks. (68) 
s 

This value is obviously much greater than the proposed lower bound given by 
Eq. (EH). 

So far, we assumed a radiation dominated universe. If the state equation is 
instead taken as p = (7 — l)pc 2 , we obtain analogously 

^ = ^Ks, (69) 

s 7 

agreeing with Eq. (|68[) when 7 = 4/3. 

Other cases analyzed in jllj . taken under widely different physical conditions, 
turned out to give values of r//s of roughly the same order of magnitude as above. 
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